SOLUTIONS OF BIGRADED TODA HIERARCHY 



CHUANZHONG LI 



Department of Mathematics, USTC, Hefei, 230026, Anhui, P. R. China 
Department of Mathematics, NBU, Ningho, 315211, Zhejiang, P. R. China 

Abstract. The {N, Af )-bigraded Toda hierarchy is an extension of the original Toda lattice 
hierarchy. The pair of numbers (N, M) represents the band structure of the Lax matrix which 
has N upper and M lower diagonals, and the original one is referred to as the (1, l)-bigraded 
Toda hierarchy. Because of this band structure, one can introduce AI + N —1 commuting flows 
which give a parametrization of a small phase space for a topological field theory. 

In this paper, we first show that there exists a natural symmetry between the (iV, M)- 
and (M, A^)-bigraded Toda hierarchies. We then derive the Hirota bilinear form for those 
commuting flows, which consists of two-dimensional Toda hierarchy, the discrete KP hierarchy 
and its Backlund transformations. We also discuss the solution structure of the (TV, M)-bigraded 
Toda equation in terms of the moment matrix defined via the wave operators associated with 
the Lax operator, and construct some of the explicit solutions. In particular, we give the 
rational solutions which are expressed by the products of the Schur polynomials corresponding 
to non-rectangular Young diagrams. 
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1. Introduction 



The (A^, M)-bigraded Toda hierarchy, denoted by (A^, M)-BTH, is an integrable system (see 
e.g. [I1I2]), and its Lax operator is given by 

£ := + ujv-iA''-^ + --- + U0 + --- + m_mA-*^. 

where N, M > 1, A is a shift operator which can be expressed as an infinite matrix in the form, 
A = {Ei^i+i)i(zz. In terms of an infinite size matrix, the Lax operator C has the band structure 
with A^ upper and M lower nonzero diagonals. The (A^, M)-BTH is then defined by 

^^^^ a£ _ I [(£"+i-^)+,£], if 7 = « = l,2,...,Ar, 

•^^7," [ [-(£"+i+^)_,£], if 7 = /3 = -M+1,...,-1,0, 

Here we call the A^ + M — 1 numbers of the flows for n = the primaries of the BTH which 
describe the small phase space of a topological field theory (TFT), and the flows with n > 
correspond to the gravitational descendants in this TFT. 

In the case of A^ = M = 1, the (1, 1)-BTH is the original Toda lattice hierarchy, and the 
primary flow is just the Toda lattice equation [3lll]. The (A^, M)-BTII with A^ > 1 or M > 1 
is then considered as an extension of the original Toda lattice hierarchy. One should also note 
that the case with infinite A^ and M corresponds to the two-dimensional Toda hierarchy where 
we have two independent Lax operators (one defined near infinity and the other defined near 
zero in the spectral space, more precisely, one considers two cases with (l,oo) and (oo,l)). 
Then the (A^, M)-BTH can be naturally considered as a reduction of the two-dimensional Toda 
hierarchy by imposing an algebraic relation to those two Lax operators (see [SHTj). 

In [2] , we showed the integrability of an extended version of ( A^, M)-BTH by writing the 
hierarchy as a bilinear identity, and introduced the r-functions. Here the extension implies 
that the hierarchy has additional logarithmic flows, and this version is called the extended 
BTH (see [HIH]). In this paper, we are interested in constructing several explicit solutions of 
the BTH, and as a first step, we only consider the non-extended version of the BTH based on 
our previous study [2]. 

The paper is organized as follows. In section [2l we give a brief summary of the original 
Toda lattice hierarchy whose Lax operator is given by a tri-diagonal matrix. We also discuss 
briefly (2, 1)-BTH as an extension of the Toda hierarchy and describe the t2fi-G.ow defined by 
the square root of the Lax operator. In particular, we mention that there exists nonlocal terms 
in the equation, and make a remark that the flow also appears in the recent paper [9]. In 
section [3], we give the explicit form of (A^, M)-BTH and the r-functions. This section is a 
brief summary of our previous paper p| without the logarithmic flows. Here we express the 
coefficient functions in the Lax operator in terms of the r-functions in the similar manner 
discussed in [T0l[T3]. We also discuss some details of the (2,2)-BTH as an example. In section 
m we show the equivalence between the (A^, M)- and (M, A^)-BTHs by using the Hirota bilinear 
equations found in [2] and gauge transformation in [Hj. To be illustrative, we also give some 
simplest concrete examples of equivalence between flows of (2, 1)-BTH and (1,2)-BTH. This 
equivalence is explicitly shown in the examples given in section O In section El we derive 
the Hirota bilinear equations for the primaries of the (A^, M)-BTH. Then we construct the 
r-functions in terms of the moment matrix defined naturally via the wave operators introduced 
in section [3] (also see [2|). In section [6l we construct rational solutions based on the r- function 
formulas derived in the previous section. In particular, those rational solutions are given by the 
products of two Schur polynomials depending on two different sets of flow parameters ta,n and 
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t[5,n in (II-IP - Contrary to the case of the original Toda hierarchy where the rational solutions 
are given by the Schur polynomials of rectangular Young diagrams, the rational solutions of the 
BTH are parametrized by non-rectangular Young diagrams. Finally, in section [TJ we summarize 
the results and give some discussions. 

2. Tridigonal Toda lattice hierarchy and generalization 

Here we briefly explain the BTH as an extension of the original Toda equation, and present 
some connection to the recent study in [H]. The main point is to explain the structure of an 
additional symmetry generated by a fractional power of the Lax matrix. The Toda lattice 
equation is written in the form. 



(2.2) ^ - n=1.2.... 

Eq.f i2.2p has the Lax representation with a tridiagonal semi-infinite matrix L by 

Bi = [L]>o , 

where [i^]>o is the upper triangular part of the matrix L given by 



(2.3) ^^-[B.Li 



(2.4) L 



/&i 1 

a2 62 1 

as 63 1 

a4 64 

V; ; ; ; 



If we consider the bounded Toda lattice equation, the Lax matrix will have finite size. 

For semi-infinite Toda equation there exists a sequence of r-functions {t„ : n > 0} with 
To = 1 defined by the a„, 6„. by the formulas 

(2.5) ttn = ^"^'""^ , bn = ^\og 



n-1 



dti 



Tn-l 



Then we can write the Toda lattice equation in the Hirota bilinear form, 
(2.6) DlTn ■ Tn = 2r„+ir„_i, 

where Di is the usual Hirota derivative. For the A;-th flow-parameter of the Toda lattice 
hierarchy, D^. is defined by 



ik—i'i- 



(2.8) 



The hierarchy of the Toda lattice is defined by 
dL 



[Bk, L], Bk = [L^] 



>0) 



k = 1,2,3,... 



The r-functions of the Toda lattice hierarchy obey the following equations 
(2.9) [Dk - Pk{D)]Tn+i ■ T„ = 0, A; = 2, 3, 4, 



where Schur polynomial Pk{D) is defined by 

oo 111 

(2.10) e^^-^^""^'" =J2Pk{D)^\ D = {D„-D„-D,,-D„...). 

When k = 2, Hirota equation becomes 

(2.11) {D2-Dl)T^+,-Tn = 0. 

Eg. (12. 11 1) together with eg. (12. 6 p will give the nonlinear Schrodinger eguation which can be seen 
as the second member of the Toda lattice hierarchy pLLj. 

A natural guestion will be how about generalized band structure of the Toda Lattice hierarchy 
which is just (1,1) tridiagonal band matrix. 

For example, for (2, 1) Heissenberg band structure of Lax matrix 



(2.12) L 



a2 





Cl 

as 




1 




1 

&4 



(2.14) 



the eguations for Toda flow, i.e. 
(2.13) dt,,,L=[L>o,L], 
will lead to the following Blaszak-Marciniak lattice eguation 

^t\X)bn CfiCLn+l 0,nCn—l 

dtxfiO^n = 0,n{bn — f'n-l)- 

Eg. (12. 141) is also eguivalent to the Bogoyavlensky-Narita lattice given in [15]. In [H], Blaszak 
and Marciniak considered the local flows which correspond to the integer powers of Lax opera- 
tors. In this paper, we construct nonlocal flows using fractional power of Lax operator. Because 
of the (2, l)-band structure, one can define the sguare root of Lax matrix which will be shown 
in detail in the next section. Using the operator La, we give a new flow 

(2.15) 

which further leads to 

- 6„ + c„(l - A)(l + A)-ic„ 
a„(l-A-i)(l + A)^ic,. 



'^t2fiL — [L^Q, L], 



(2.16) 



dt2flbn 



After denoting V, as eg.( l2.16p can be rewritten as 



(2.17) 



After transformation 



9t2,obn 
^t2.o^n 



bn+1 



Ojrt 



(2.18) 
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eq. (12.16p becomes 

{dt2,o'^n+l + dt2^o^n = bn+1 — + c\ — C^j^^ 
dt2fiK = fln+l - On 

which is just eq.(40) in fT5] and also related to the system (10)-(12) proposed in [T6j. In [T5] . 
Svinin gives general constructions of such flows. Studying the relation between two ways of 
constructing nonlocal flows, i.e. the way in [15] and the way in this paper might be an interesting 
problem. 

Eq. ( I2.14P and eq. (I2.17P are in fact ti flow and t2 flow of the case n = 3, a = — 1 without 
center extension in [9j, i.e. solutions do not depend on y variable. In this paper, we directly 
introduce a new fractional Lax matrix instead of using Casimir construction used in P for 
constructing the Lax equations. We also generalize these results to (A^, M)-band matrix. For a 
flnite-sized Lax matrix, its fractional power may not be well-deflned. This leads to a difficulty 
to give symmetric flows generated by fraction powers of Lax matrix. But for a bi-inflnite band 
matrix, one can deflne the fraction powers and further deflne other additional flows which 
commute with the original Toda flow. This generalization leads to the BTH which might be 
also seen as a general reduction of the two-dimensional Toda lattice hierarchy. In the next 
section, we give the continuous interpolated version of BTH [2] and later present the matrix 
version of the BTH in bi-inflnite and semi-inflnite band matrices. 

3. The bigraded Toda hierarchy (BTH) 
The Lax operator of the BTH is given by the Laurent polynomial of A pLj 

(3.20) £ := + UN-iK^'^ + ■ ■ ■ + + ■ ■ ■ + M-a/A"*^ 

where N,M > 1, A represents the shift operator with A := e'^^"' and "e" is called the string 
coupling constant, i.e. for any function f{x) 

Af{x) = f{x + e). 

The C can be written in two different ways by dressing the shift operator 

(3.21) C = VlA^Vl^ = VrA-^^V^^, 
where the dressing operators have the form, 

(3.22) Vl = 1 + wiA~^ + W2A-'^ + ..., 

(3.23) Vb^ = Wo + WiA + W2A^ + .... 

Eq. (13.211) are quite important because it gives the reduction condition from the two-dimensional 
Toda lattice hierarchy. The pair is unique up to multiplying Vl and Vr from the right by 
operators in the form 1 + aiA~^ + a2A~^ + ... and clq + diA + 02 + . . . respectively with 
coefficients independent of x. Given any difference operator A = J^k^^A'', the positive and 
negative projections are defined by = ^fc>o^fcA'^ and A„ = J2k<o^kA''. 

To write out explicitly the Lax equations of BTH, fractional powers and £s" are defined 

by 

C^=A + J2 ^kA\ =Y. 

fc<0 fc>-l 
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with the relations 

{C^rf = (£ir)*^ = C. 

Acting on free function, these two fraction powers can be seen as two different locally expan- 
sions around zero and infinity respectively. It was stressed that £jv and jC,m are two different 
operators even if = M{N, M > 2) in [1] due to two different dressing operators. They can 
also be expressed as following 

Let us now define the following operators for the generators of the BTH flows, 

if 7 = a = l,2,...,A^ 
if 7 = /3 = -M + 1,...,-1,0, 



(3.24) B. 



7,n 



Definition 3.1. Bigraded Toda hierarchy (BTH) in the Lax representation is given by the set 
of infinite number of flows defined by 

, , dC, \ [(5a,n)+,/:], if 7 = a = l,2,...,iV, 

3.25 = <^ 

\\-{Bfi^n)-.C.l if 7 = /3 = -M + l,...,-l,0. 

We need to remark that this kind of definition is equivalent to the definition in [1] which is just 
a scalar transformation about time variables using gamma function. The original tridiagonal 
Toda (i.e. Kostant-Toda) hierarchy corresponds to the case with N = M = \. 

One can show [2] that the BTH in the Lax representation can be written in the equations of 
the dressing operators (i.e. the Sato equations): 

Theorem 3.2. The operator C in fl3.20p is a solution to the BTH fl3.25p if and only if there 
is a pair of dressing operators Vl and Vr which satisfy the Sato equations, 

(3.26) d^,nVL = -{B^,n)-VL, d^,nVR = {B^,n)+VR 

for -M + 1 <'j < N andn>0. 

The dressing operators satisfying Sato equations will be called wave operators. By wave 
operators we will give the definition of tau function for BTH as following. 

According to paper p], a function r depending only on the dynamical variables t and e is 
called the tau-function of BTH if it provides symbols related to wave operators as following, 

(3.27) P.: = l + !^ + ^ + ,,,:^ ^<^-'-|ff-^) . 

A T{x,t;e) 

-llN. 



p-1. ^ 1 ._ T{x + e,t+[X-Y;e) 

^ ■ A A2 ■ T(x + e,t;e) 



(3.28) 

(3.29) Pr: = Wo + w^X + W2X^ + 



T{x + e,t+[Xf';e) 



t(x, t: e 



5 ^5 

M. 



(3.30) P^^: = w'o + w[X + w',X^ + ...:= ^^'';^ \'\ 



where [A and [A]^^ are defined by 



^ ^-'^ ^0, 7 = 0,-l----(M-l), 



7 = iV,iV-l,...l, 




For a given pair of wave operators tlie tau-function is unique up to a non-vanisliing function 
factor which is independent of x and t^^n with all n > and — M + 1<7<A^. 
Then we get 



(3 31) P - V Pn{-dL)r{x,t] e) p-i Pn{dL)T{x + e,t]e) 



(3.32) P« .- A , Pn -1^ r(x + 6,t;e) ^ ' 

where P„ are the elementary Schur polynomial as defined in (12. 101) . Here the operators di and 
Br are defined by 

dL=[ — ^ rrdt ■.l<a< N 

Or = \ ^dt, : -M + 1< /3 < 

The dressing operators Vl and Vr can be expressed by function T{x,t; e): 

An T-)-! \^ \nPn{ dR)T{x,t;e) 

(3.34) A , -)_^A ^(^+ , ) ■ 

One can then find the explicit form of the coefficients Ui{x,t) of the operator C in terms of 
the r-function using eg. (13.21 p as 



(3 35) u {x t) = P^-i(P'L)r{x + (z + l)e, t; e) o t{x, t; e) ^ Pm+;{Dr)t{x + e, t; e) o t{x + ze, t; e) 

r(x, t; e) r(x + (z + l)e, t; e) r(a:, t; e) t(x + (z + l)e, t; e) 

where and Dr are just the Hirota derivatives corresponding to di and Sr respectively. 

The BTH can be also written in the matrix form with the identification of the shift operator 
A as the infinite matrix having zero entries except I's in the upper diagonal elements and all 
other functions about x as diagonal infinite matrix [6] . But now we only consider its reduction, 
i.e. the corresponding semi-infinite matrix form in the following. Then we rewrite the coefficient 
Ui{x,t) as Uij{t) and rewrite t{x + e,t) as Tj{t). We can find the corresponding semi-infinite 

matrix forms Vli'Pl^ ,'Pri'Pr^ corresponding to Vl.VJ^^ ^Vr,V^^ respectively as following 
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(3.36) 



Vl 



1 














1 





n 






P2(-9l)t2 


PA 


-dL)T2 


1 










P3(-9l)t3 




-9i)T3 


Pi{-dL) 


T3 




T-3 


T3 



/ 1 

Pi(8z.)ri 



(3.37) 





1 






1 



P2{dL)Tl Pl{dL)T2 
Pz{dL)Tl P20L)r2 Pl{dL)T2 







1 







(3.38) 



(3.39) 



V 



V 


Tl 


T2 


— ■ — 1 

T2 












J 


21 


Pi0r)ti 


P2(dR)Tl 


ft(9fl)Tl 






TO 


TO 


TO 


TO 
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Pl{dR)T2 


P2(9fl)T2 








Tl 


Tl 














a 


Pi{dR)r3 






T2 


T2 















a 












T3 














• y' 
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Pl{~dR)Tl 


P2(~dR)T2 


-P3(-9j{)-n3 




Tl 


T2 


T3 


T4 







H 


Pi(-9h)t^ 


P2(-aj{)T3 






T2 


T3 












Ik 


Pi(-9j{)-n3 




T3 


T4 













23. 










T4 





After the following transformation ttj ^ = O'jj+i, the matrix representation of C can be expressed 
by {ai,j)i,j>i with 



(3.40) a,,,(t) = 



Pi-j+N{DL)Tj O Ti-i _ Pj^i+MiDR)Ti O Tj_i 



Note here that those expressions immediately imply 

ttij = 0, if j < -M + i or j > N + i. 

That is, the Lax matrix C has the (TV, M) band structure. 

As an example, we will give some concrete results on (2, 2)-BTH in the following subsection 
from which we can see some general patten of {N, M)-BTH. 

3.1. Example of the (2,2)-BTH. Let us summarize this section by taking the (2,2)-BTH. 

The Lax operator is 

(3.41) L = A^ + uiA + uo + u_iA"^ + «_2A"^ 

i_ 1 

Then there will be two different fraction power of L, denoted as and respectively as 
following form 

(3.42) = A + ao + a^-^A'^ + a_2A"^ + • • • , 



(3.43) Lfj = a'_^A-^ + + a[A + a'^A^ + .... 

We can get some relations of {af, i < 0}, {a^-; j > —1} with {uf, —M < i < N — 1} as following 



(3.44) 



ao(x) = {l + A)-^ui{x), 



"-1 



,(H-A-i)-ilog«-2(a:) 



Then by Lax equation, we get the ^2,0 flow of (2,2)-BTH 



(3.45) 

which correspond to 



(3.46) 



^2,0^^0 (a;) 
^2,0^^-1 (2:) 

^d2,0U-2{x) 



d2,oL= [A + {l + A)'\^{x),L] 

uo{x + e) — ^0(3;) + ui{x){l — A)(l + A)^^ui{x) 
u_i{x + e) — u_i{x) 

u_2{x + e) - u^2ix) + ii_i(x)(l - A-^)(l + A)-^ui{x) 
u^2{x){l - A-^){1 + A)-'ui{x). 



From eqs. 03.461) . we can find the equations have infinite terms because of (1 + A)~^ which comes 
from the fraction power of the Lax operator. Just like the method in [T7j , for avoiding infinite 
sums, we use auxiliary function aQ{x) with which we can rewrite eq. fl3.46p as 



(3.47) < 



'^2,0^^1 (a:) 
diflUoix) 
d2,ou^i{x) 

d2,0U-2{x) 



Uo{x + e) - uo{x) + ui{x){ao{x) - ao{x + e)) 
u^i{x + e) — u^i{x) 

u-2{x + e) - ^-2(2^) + M-i(a;)(ao(x) - ao(a; - e)) 
U-2{,x){aQ{x) — ao(x — 2e)) 



(3.49) 



,d2flaQ{x + e) + d2flaQ{x) = uq{x + e) - uq{x) + ui{x){aQ{x) - ao(x + e)). 

The tifl flow will have finite terms as following because it does not use the fraction power of 
Lax operator L, 

(3.48) diflL = [A' + uiA + mq, L] 

which correspond to 

' diflUi{x) = u_i{x + 2e) — u^i{x) 

difiUQ^x) = u_2{x + 2e) — u_2{x) + Ui{x)u^i{x + e) — u„i{x)ui{x — e) 
diflU^i{x) = Ui{x)u^2{x + e) — M_2(a;)tti(x — 2e) + u^i{x){uq{x) — uq{x — e)) 
^difiU^2{x) = M_2(x)(uo(x) - uo{x - 2e)). 

For t_i flow, equations will also be complicated because of another fraction power of L. The 
equation is 

(3.50) d^iflL = -[e(i+'^"')"''°g"-2A-\L] 

which corresponds to 

'5-1,0^^1 (2:) 

d-iflUo{x) 

5-i,oW-i(a;) = e(i+^^')"'i°s"-2(^')(uo(x) - uo{x - e)) 

^a_i,oU-2(a;) = u_i(x)e(i+^"')"'i°s«-2(x-.) _ g(i+A-i)-iiog«_2(x)^_^^^ _ 



g(l+A-l)-llog«_2(a;+2e) _ g(l+A-i )"! log «_2 (x) 



(3.51) < 



ui[x)e 



(l+A-i)- 



;m_2(x) 



Ui[x — e) 



From eqs. dS.Sl]) . we can find the equations have finite terms but every term has infinite multi- 
phcation because of factor e^^"*"^ ^ iog«_2(x) -y^^j^jcj^ comes from the square root of Lax operator. 
Similarly for avoiding infinite multiphcation, we use auxiliary function a'_^{x) with which we 
can rewrite eq. fl3.5ip as 

d^iflUi{x) = a'_i{x + 2e) — a'_i(x) 

d-i,oUo{x) = ui{x)a'_-j^{x + e) — a'_i{x)ui{x — e) 

(3-52) )d-iflU-i{x) = a'„i(2;)(Mo(x) - wo(a; - e)) 

d-i^oU-2{x) = U-i{x)a'_i{x — e) — a'_i{x)u-i{x — e) 
9_i,oa.'_i(a^) = U-i{x). 

Infinite sums or infinite muhiphcations are important properties of BTH because of nonlocal 
operators. For finite Lax matrix, it is not easy to construct its fraction power. That is why we 
do not use fraction power of finite matrix to give Lax equations. 

4. Equivalence between (N, M)-BTH and (M, A^)-BTH 

In this section, we prove that there is an equivalence between (A^, M)-BTH and (M, A^)- 
BTH in three ways. One is to prove the equivalence in the Hirota bilinear identities basing 
on [2]. The second one is for equivalence in specific Hirota bilinear equations. At last, we will 
prove the equivalence between their Lax equations using transformation in fT?]. To see the 
equivalence clearly, one explicit example, i.e. equivalence between (1,2)-BTH and (2, 1)-BTH 
in Lax equations under transformation will be shown in detail. 

4.1. Equivalence in the Hirota bilinear identities. Firstly after denoting r{x,t) as t{x — 
we get the following Hirota bilinear identity p], i.e. for each m G Z, r G N, 

ResA {A^^+'"-V(x,t - [A-^) X r{x - (m - l)e,t' + [A-i]^)e«^(^'*-*') 

(4.1) = ResA {A-""^+™-V(a; + e,t+ [Xf ) x t{x - me,t' - [A]^Oe-««(""''*-*') 

where 

N 

^d\t) = $:$:A^("+-^)f.,., 

ra>0 a=l 
n>0 /3=-A/+l 

The most important property that BTH has is that there are Nr and Mr in both sides of 
HBEs (14. ip . These two terms show the principal difference of BTH from the two-dimensional 
Toda hierarchy, i.e. the constraint (I3.2ip . 

Hirota bilinear identity eq. (l4.ip can lead to the following identity under the transformation 

m I— — m, X X — me, A i— ?► A^^ 

ResA I A-U-^^+™r(x - me, t - [A]^) x r{x + e,t' + [Xf)e^^^^-'~''^ 



(4.2) = ResA |A-iA*^^'+'"r(x - (m - l)e, t + [A" x ^i^, t' - [A-i]^^)e-««(^"''*-*') 
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where 

N 

n>0 a=l 
n>0 /3=-A/+l 

The eq. (l4.2p can be rewritten as following identity after the interchanging of t and t' 
ResA |a-^A-^"+"^t(x + e,t+ [Xf) x t{x - me,t' - [A]^)e«^(^'*'-*)} 

(4.3) = ResA |a-U*^"+"^t(x, t - [X-^] x r(x - (m - l)e, t' + [A^YOe"^^^^"''*'"*^ 
which can be further rewritten as following 

ResA \x-^X^'^'^+"'T{x,t- [X-^) X ^i^ - (m - l)e,t' + [A-i]^^)e««(^'*-*') 



(4.4) = ResA { A-U~^^'+'"t(x + e, t + [A]^) x r(x - me, t' - [Xf)e~^^^^'"^'-''^ } . 

Eq. (l4.4p is obviously (M, A^)-BTH comparing to eg. (14. II) if we change time variable t^^„ to 
ti-^^ni i-e. subscript L ^ R. Therefore for (M, A^)-BTH, Eq. (l4.4p is in fact changed into the 
following equation under transformation „ — )■ ti_^,ri 

ResA |A-U*^^'+'"r(x,t - [A-^') x r(x - (m - l)e,t' + [X-T)e^'^'''~''^} 

(4.5) = ResA |A-^A-^'^+'"r(x + e,t+ [Xf) x t{x - me, t' - [A]^)e-«^(^"''*-*')} . 

Because transforms m i— )■ — m, x i— )■ x—me, X ^ X~^ do not change the equation itself. Therefore 
we can say that the {N, M)-BTH is equivalent to the (M, A^)-BTH under the transformation 

4.2. Equivalence in the Hirota bilinear equations. The bilinear identity eq. fl4.ip of the 

BTH can be equivalently expressed as [2j 

ResA {A^''+--V,_(^_i)(t + y + [X-^) r,{t - y - [A-^) e«-(^'-2^)} 

(4.6) = ResA {A-^^'-+™-V,+i(t - y + [A]^^) T,_„(t + y - [Xf) e^-^^'^^^y^] . 

To be specific, we will give the equivalence from the concrete Hirota equations between 
[N, M)-BTH and (M, A^)-BTH as following which will also be used to derive specific primary 
Hirota equations of BTH in the next section. 

For (iV, M)-BTH, in term with H vlli^vlli, ■ ■ ■ V^is U vLv/U ' ' ' "^+1 < A < 0, 1 < 
ttj < ^ ■, the Hirota equation is as 

v=l i[=0 i'p=Oq=l 1' ^9 1'' 

P D^'^ ^1 ks s ,^ ( n\k„-i„ 

nZJiAfSr \^rj [^^^iMIlI^I^p , (f) 

v=l n=0 is=Oq=l ^ 'i'' 
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(4.7) 

where Pk are Schur polynomial as defined in f l2.10p . After the transformation m t— )■ —m,n t— >■ 
n — m, the identity eq. (l4.7p becomes 

v=l ji=0 is=Oq=l ^ ^ 

^ ( n \>'^ '''^ ^ (-Dr ,>Y'i 9*=i-*', 

(4.8)r„^m+ir„. 

After doing the transformation D^ i = D[_^ i, —M + 1 < j < N, eq. fl4.8p becomes 

P (]-)> \K ^1 ■'I'D' Yi ( 0\kq-i„ 

v=l ii=0 i^=Og=l ^ 'S^' 



V l—av,lv' 



11=1 



(4-9) ( ■ ■ ■ 5Z n (U _ .vy^E^^i _i/ )+Mr+m(^L))^n-™+l^n- 



which is the term with ]\ y'l^^ . . . y%^^,^ [1 • • • Vi-o^as^^ < 1 - A < 
M, -A^ + 1 < 1 - < for (M, A)-BTH. So there is a correspondence of (A^, M)-BTH and 
(M , A^)-BTH in Hirota equations under the meaning of following derivatives' correspondence 
O -M + 1 < 7 < A, i.e. D'l ^ Dr, D'^ ^ Dl. 

Because (A^, M)-BTH and (M, A^)-BTH are equivalent, we can only consider BTH in the 
case of < M later. 

4.3. Equivalence in the Lax equations. Using the gauge transformation and linear transfor- 
mation mentioned in [l^, we can prove the equivalence between (A^, M)-BTH and (M, A^)-BTH 
under the meaning of Lax equations. 

Firstly we need to introduce the following proposition basing on Theorem 6 in 



Proposition 4.1. If C satisfies Lax equations f l3.25p and let $ = $('u,t) satisfies condition 

(4.10) ai^_„<l> = (S^,„)o$, -M + l<7< Ar,n>0, 

(where subscript denotes the projection to term of hP), then C = ^^^C^ satisfies the hierarchy 

(4.11) a^,„£= [(5^,„)>i,£], -M + l<7< A,n>0, 
where 

(4.12) = 

Proof. Because of f l3.25p and f l4.10p . the following calculation holds 
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= 0. 

Then we finished the proof of this proposition. □ 

Proposition 14. ll tells us that the gauge transformation from Theorem 6 in [T^ can be extended 
on the fractional powers of Lax operators. Now let us introduce the following notation for anti- 
involution map "t": If £ is as form fl3.20p . then 

C'^: = A*VM(x)+A*'-VAm(x) + --- + A-^+V-i(a:)+A-'^ 

= u_m{x + Me)A^^ + u_M+i{x + (M + l)e)A^'^-i + ■ ■ ■ + un-i{x - {N - l)e)A-^+i + A^^. 

By calculation, one can prove the following two lemmas similar as [T^ . 

Lemma 4.2. For —M + I < < N,n > 0, following identity holds for any integer k 

(4.13) {{B,,n)>ky = ((S^,„)t)<-fc. 

Using above lemma, we can prove the following lemma directly. 

Lemma 4.3. dt^ „C = [(S^,„)>o,£] can leadtod^ ^C^ = [{{B^^nV)>i, C^] andd^ ,^ = [(S^,„)>i,£] 
can lead to dt^ ^^C^ = [{{B^^nV)>o, C^. 



By above two lemmas and Proposition 14. H we can prove the following important theorem. 
Theorem 4.4. Lax equation fl3.25l) of {N, M)-BTH with Lax operator 

(4.14) Ln,m = A^ + UN-iA^'^ + --- + U0 + --- + U-mA'^ 
is equivalent to {M, N)-BTH with Lax operator 

(4.15) Cm,n = A^-^ + ua/-iA^^-i + ... + U0 + --- + U-nA'^ 
under Miura map: Uj{x,t) = U-j(x + je)e'>^-^~^' 

Proof For Lax operator Cn,m = A^ + un^iA^'^ + ■■■ + % + ■■■ + M-a/A"*^ of (A^, M)-BTH 
which satisfies Lax equations (13. 25 p . one can choose \E' = \i/(M,t) = e'-^"'^ "-mC^^,*) which 
satisfies condition 

(4.16) 5t^,„^ = (5^,„)o^, -M + l<7<iV,n>0. 

Then Cn,m = ^"^^at.a/^ = unA^ + un^iA^'^ H h wo H h A"^"^ satisfies the hierarchy 

(4.17) dt^^^CN,M = [{B^,n)>l, Cn,m], -M + 1 < 7 < iV,n > 0, 
where 

(4.18) B^^n = ^"^-B^,n^, Ui{x) = ^'^{x)ui{x)'i'{x + ie). 



Using Lemma [4. 3 [ one can derive 

(4.19) dt.,J^j,^M = [iiB,J)>o,C^^^M], -M + 1 < 7 < iV,n > 0, 

One can choose Lax operator Cm,n = -^jy = A^'^ + ■uai-iA^"'^"^ + ■ ■ ■ + Uq + ■ ■ ■ + W-atA"^ of 
(M, A^)-BTH with 

Uj{x) = u^j{x + je) = "^'^{x + ie)u_j{x + je)\E'(x) = -u_j(x + je)ei-A-Af "-^/(i^.t)^ 
The (M, A^)-BTH with this Lax matrix will be equivalent to the original (A^, M)-BTH. □ 

13 



To be illustrative, we will give some simplest concrete examples of equivalence in the next 
subsection which includes nonlocal flows of (2, 1)-BTH and (1,2)-BTH. Here we will also con- 
sider the interpolated form of BTH. At this time, we will use a powerful tool called gauge 
transformation to prove that equivalence. 



4.4. Equivalence between (1,2)-BTH and (2, 1)-BTH. Using Proposition 14. ![ we will see 

the equivalence between (1,2)-BTH and (2, 1)-BTH in detail. Here we only give the primary 
flows of them to see the equivalence. 

(1, 2)-BTH: The Lax operator of (1, 2)-BTH is as following 

(4.20) £i,2 = A + tio + M-iA"^+M_2A-2. 
(1, 2)-BTH has the following primary equations 

(4.21) ai,oA,2 = [A + %,A,2], 
and 

(4.22) 9_i,o£i,2 = -[e(i+^-')-^^°^"-A-\£i,2], 
which further lead to 

(4.23) I 9_i,oM-i(x) = e(i+'^"')"''°g"-2(^)(Mo(x) - uo{x - e)), 

[a_i,oM-2(x) = M_i(x)e(i+^'')~'i°g"-2(^-^) - e(i+^"')"'i°g"-2WM_i(x - e). 

and 

{difiUo{x) = ■u„i(x + e) - 'U_i(x), 
diflU-i{x) = U-2{x + e) - u^2{x) + u^i{x){uq{x) - Uo{x - e)), 
diflU-2{x) = u^2{x){uo{x) - Uo{x - 2e)). 

(2, 1)-BTH: The Lax operator of (2, 1)-BTH is as following 

(4.25) £2,1 = A^ + WiA + Mo + M_iA~\ 
The equations fl3.25p in this case are as follows 

(4.26) 92,o/:2,i = [A + (1 + A)-ini(x), £2,1], 

(4.27) 81,0^2,1 = [A^ + UiA + uo, £2,1] , 
which further lead to the following concrete equations 

{d2flUi{x) = ui{x + e) - ui{x) + ui{x){l - A)(l + A)~^-ui(x), 
d2,oUo{x) = ■u_i(x + e) - u-i{x), 
d2,oU-iix) = u_,{x)il - A-i)(l + A)-'u,{x), 

{difiUi{x) = -u_i(x + 2e) - n_i(x), 
difiUo{x) = u^2{x + 2e) — u^2{x) + Ui{x)u^i{x + e) — U-i{x)ui{x — e), 
difiU-i{x) = u_i{x){uo{x) - Uo{x - e)). 

It seems that eq. fl4.28p and eq. fl4.29p are quite different from eq. (14.231) and eq. (14.240 respectively. 
In fact after doing gauge transformation on (2, 1)-BTH as following, we can find the equivalent 
relation between these equations. 
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Now we consider function has form 

(4.30) ^_g(i-A)-iiogu_i(x)^ 

then £2,1 := 0^^£2,i0 will have form 

(4.31) £2,1 = ^^2A^ + ^lA + ^^0 + A"^ 
The relation of Vi{0 < i < 2) and Ui{—1 < i < 1) are like 

(4.32) V2 = 0"V(a; + 2e), t;i = 0~-^mi0(x + e), vq = uq. 
Therefore we get following new flows on new Lax operator £2,1 using Proposition 14.11 

(4.33) 92,o£2,i = [e(i+^)"'i°s^^(^)A,£2,i], 

(4.34) ai,o£2,i = hA' + t^iA\£2,i], 
which further leads to 

r ^2,0^^0(2:) = e(l+^)"''°S^2{a=) _ g(l+A)-Mog^2(x-e)^ 

(4.35) I d2,oVi{x) = e(i+^)"'i°s''2(^)(i;o(a; + e) - vo{x)), 

[•92,0^^2(2:) = vi{x + e)e(^+^)"'^°s^2(x) _ g(i+A)-iiogi,2(x+E)^^(^^)_ 

and 

(diflVoix) = vi{x) -vi{x-e), 
diflVi{x) = V2ix) -V2ix - e) + vi{x){vo{x + e) - Vo{x)), 
di,oV2{x) = V2{x){vo{x + 2e) - Vo{x)). 

Comparing eq. (l4.35p . eq. fl4.36p with eg. (14.231) and eq. fl4.24p . we can find these two pairs of flows 
are equivalent under Miura map Uj = v_j{x + ie)(0 < j < 2). This proves the equivalence 
between (1,2)-BTH and (2, 1)-BTH. 

Using bilinear identities got by comparing every term of Hirota bilinear identities of BTH in 
this section, in the next section we will derive all primary Hirota equations of BTH to see its 
inner structure. 

5. Hirota equations and solutions of the BTH 

Hirota bilinear equations(HBEs) are central object in Sato theory. From HBEs, we can 
derive the structure of solution. This is a great motivation for us to consider HBEs of the 
BTH. The Hirota bilinear equations of the BTH can be derived from (14.60 which comes from 
HBEs in [2j. In particular, the following proposition will list all the Hirota equations for the 
primary variables, i.e. t^^„ with n = 0. 

Proposition 5.1. For {N, M)-BTH, we have the following identities for primary derivatives 
which are equivalent to all the primary Hirota equations. 

(5.1) (d^,o - PmMDr)) Tn+i o r„ = 0, 

(5.2) (-D/3,o^-A/+i,o - 2PM+/3+i(^i?)) T-„ o r„ = 0, 

(5.3) D/s^oDNfiTn o r„ = 2PM+,3-l(^i?)T:«+l o 

(5.4) -Da,Op'-A/+l,0^n ° = '^Pn -a{D L)Tn+l O T„_i 
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(5.5) (^DaflDN,0 - 2PN-a+2{DL)^ T„ O Tn = 0, 

(5.6) (p^fl - PN-a+l{DL)) r-a+l O T„ = 0. 

where are Schur polynomial as defined in fl2.10p . 



Proof. In eq. fl4.6p . for term y^^, —M + 1 < /3 < in (A^, M)-BTH, following Hirota equation 
holds 



1=0 



(5.7) 

which can also be got from eq. fl4.7p . 

For 1 < a < in (A^, M)-BTH, following Hirota equation holds 



i=0 



(5.8) 

For 7/^,0, + 1 < /3 < in (A^, M)-BTH, the following Hirota equation holds 

^ (— -D^i)* 2}~^ 

which is 



(2-^Af(l+^)+Mr-m(-^-R) ~ D ii^QPMr~m{D p!))Tn+lTn-m — D p fiPMr+m{D l) '^n—m+l'^n- 

Set r = 0, m = —1, for term with y/j^o, —M + 1 < /3 < 0, we get 

(2-PAf+^+l(^/?) — DiiflD_M+lfl)Tn+lT~n+l = 0. 

When r = 0, m = 0, 

for term with -M + 1 < /3 < in (A^, M)-BTH, we get 

(5.9) {Pm+p{Dr) - Dp^^T,,+^Tn = 0. 

Set r = 0, m = 1, for term y^,o, + 1 < /3 < in (A^, M)-BTH, we get 

2PM+li-l{D Il)Tn+lTn-l = D p fiD M fiTnTn- 

For term ^q, 0, 1 < a < A^ in (A^, M)-BTH, following Hirota equation holds 

^ /I _ ■\\^N{l-^){l~i)+Nr+mi^L))Tn~m+lTn = D afiPMr-m{D RjTn+lTn-m-, 

(5.10) 

which is 
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For r = 0,m = — 1, it further leads to 

--M+lflTn+lTn+l- 

Set r = 0, m = 0, for term with ya,o, 1 < a < in (A^, M)-BTH, following equation succeeds 

(^Af-(Q-l)(-DL) - Da,o)'''n+lTn = 

Set r = 0, m = 1, we get equation 

{-^DafiDNfi - PN-a+2{DL))rnTn = 0. 

You can get the primary Hirota equations from the value of A; = 0, 1. The other values of k 
will give higher order derivatives of the hierarchy. 

When k = 0, the equations eq. fl5.7p and eq. fl5.8l) will give 

(5.11) 

PNr+m (Dl) PMr-m{D RjTn+lTn-m- 

when r = 1, m = 0, eq. fl5.lip becomes 

(5.12) PN{DL)rn+irn = Pm{Dr) 
when r = l,m = —1, eq. fl5.1ip becomes 

(5.13) PN-iiDL) 
when r = 1 , m = 1 , we get 

(5.14) P^+i(I^L)r„r„ = Pm-i{Dr) 

The three equations eq. fl5.12p . eq. (l5.13p . eq. fl5.14p can be derived from eq. (15. ip -eq. (15.40 . These 
equations discussed above are all the primary Hirota equations. The other values of r and m 
will include higher derivatives of the hierarchy. □ 

From these primary Hirota equations, we can see BTH has ample structure information. We 
can say BTH contains discrete KP equation(from eq. (l5.2p and eq. (l5.5p ). NLS equations (from 
eq. (l5.ip and eq. (l5.6p ) and 2-dimensional Toda lattice equation(from eq. (l5.3p and eq. (l5.4p ). 
Comparing with Hirota equations of the two-dimensional Toda hierarchy we find BTH have 
more constraints on equations which comes from the equivalence of dt^ „ and dtQ „ . These 
information help us to get the solution of BTH which will be given in the next subsection. 
From all the primary Hirota equations mentioned above, we can get that the solution of BTH 
should have double- wronskian structure [T7]. In paper [T^, if we impose the vanishing of dy 
derivatives on tau functions the molecule equation in fact becomes our (2, 1)-BTH. Using the 
same method in [T7], we proved the double- wronskian solution structure satisfy all the primary 
Hirota equations of {N, M)-BTH. But later we find there is another much simpler way to get the 
structure naturally which will be mentioned in the next subsection. So we prefer this simpler 
way to the way in [17j . 
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5.1. Solutions of the BTH in the semi-infinite matrix. Now we construct tau function 
for the BTH in the semi-infinite matrix representation, that is, (oj,j)i,j>i- In order to do this, 
we first introduce the wave operators Wl and Wr associated with the dressing operators Vl 
and Vr, 



N 



(5.15) Wl{x, t, A) = Vl{x, t, A)oexpiJ2Y. A^^^+'^^^^t., 

\n>0 a=l 

(5.16) WRix,t,A)=VRix,t,A)oexp(~Yl Yl A-*^("+i+^)t;3,„ 

\ n>0 /3=-A/+l 

By sato equations f l3.26p . we can have following identities [2] 

[Ba,n)+m, a = N...l, 
-(S„J_>Vl, a = M + 1, 



(S„,„)+W,j, a = N...l, 
-(S„,„)_>Vfl, a = M+1. 

One can then prove that the product Wj^^Wr is invariant under all the fiows, i.e. 
Therefore 

Wj:'WR{x,t,A) = Wz'Wr{x,0,A) = V£'Vr{x,0,A). 
This implies 
(5.17) iVz'VRm 



N 



= E E A"^^"-^^-"^^t.,n o {VE'VRm o exp 5^ A-''^'"''^^^ t, 

\n>OQ=l / \n>0/3=-A/+l 

If we let Tq = l,Ti = 0{i E Z_); then bi-infinite matrix everywhere will become semi-infinite ma- 
trix. Representing the product V^^Vr in the semi-infinite matrix, i.e. V^^Vr and considering 
identity f l3.37p and identity fl3.38p . eq. fl5.17p can be written as following 

(5.18) {VE'VRm 



( 1 Pi(t„) P2{t^) ... \ 



1 Pi(t„) ... 
1 ... 



V 



^l'^k(O) 



/I ... \ 

Pi{tp) 1 ... 

p,{^) p,{tp) 1 ... 
V / 



where Pk are Schur polynomial as defined in f l2.10p . To construct tau functions, we define the 
moment matrix Moo(t) as, 

M^it) := {Vl'VR){t). 

Proposition 5.2. The constrained condition of BTH is equivalent to matrix M^o satisfies 
identity 

(5.19) A'^^Moo = MooA- 
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Proof. Using constraint on C, i.e. dt^^C = dt^^C, we can get dt^^M^o = dt^^M^ which can 
lead to eq.(EII2]). ' ' ' ' □ 

Direct calculation can lead to following proposition. 

Proposition 5.3. The matrix satisfies identity 

Let Mi be the i x i submatrix of Moo of the top left corner. By eq. (l5.18p . each r- function 
can be obtained by the determinant [12] , 

Ti{t) = det {Mi{t)) . 

Now, we will consider the detailed structure of tau functions of BTH from the point of reduction 
of the two-dimensional Toda hierarchy. 

As we all know, the tau functions of the two-dimensional Toda lattice hierarchy are given by 





Co,o 


Co,i 




(5.20) ri = 




^1,1 . 


C*l,i-1 






Cj-1,1 • 


• ^-1,4-1 



where 



k,l=0 

We should note here that the coefficients Cij^k,i are totally independent. 

As the original tridiagonal Toda lattice is (1, 1) reduction of the two-dimensional Toda lattice 
hierarchy. Therefore to get the solution of tridiagonal Toda lattice, we need to add factor 
S{X — /i) under the integral in the definition of Cij, i.e. 

(5.21) J J p{X,fi)6{X- fi)X'fi^e^^=o^"^"+J:u^oynf^''dXdi^, 
which can further lead to 

(5.22) J p(A,A)A'+^eS-o(-"+?/")A"rfA. 

After changing x, y time variables to ta, t^, eg. (15.221) become a new function 

p(A,A)A^+^e«^(^'*")+««(^'*'')rfA 



which corresponds to (1, 1)-BTH. 

Denote un and um as the N-th root and M-th root of unit. For {N, M)-BTH, new function 
Cjj(new form of Cij) have the following form 




p{X,fi)6{X^ - /OA>^e«^(^'*")+««(^'*'')dAd/i 
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N-1 M-1 



p=0 q=0 

where 

(5.23) tl^ 

(5.24) = «)^^("+^+^)t,,„. 

Transforms f l5.23p and fl5.24p can be seen as the twist of exponential solutions because of the 
bigraded structure. 

Because in the next we will consider rational solutions of BTH, we write Cij further into 

N-l AI-1 oo 
= X] X] X] '^iip,q^kita)Pl{tf}), 



p=0 q=0 k,l=0 



where 



a 



pKA^,a;l,Ai^)(a;^A^)'K,Ai^)^'KA^)'=(a;?,Ai^)'rfA. 



We can find coefficients {c^jp i,j,k, I > 0;0 < p < N — 1,0 < q < M — 1} satisfy 
(5.25) c. 



k+N,l 



k,l+M 



which tells us that Pm{ta)Pn+M{ii3) and Pm+N{ta)Pn{ti3) always appear at the same time. So 
the element in position of k row and / column in initial moment matrix can have the following 
form 

N-l M-1 

(5-26) (^^"'^^(O)),. = EE<oi;;'- 

Therefore tau functions of the BTH can be explicitly written in the form 





Co,o 


Co,i 




(5.27) n = 


Ci,o 










Cj-1,1 . 





With the definition of Cj j, in the next part we will construct Lax matrix solution using 
orthogonal polynomials which are nothing but wave function in semi-infinite vector form. Before 
we do that, we need some formula about property of Schur Hirota derivatives described in the 
following lemma 



Lemma 5.4. Schur derivatives have following formula 



(5.28) 



Pn{DL)T, ■ = ^ Pk{dL)n X PK-^l)^,, 



k+l=n 



(5.29) PiidL)[0,l,2,...,i-l]L 

(5.30) PK-^L)[0,l,2,...,z-l]i 



(5.31) 



Pn{DR)Ti ■ Tj 



[0,1,2, ...,^-2,Z + /-l]i, 

(-l)'[0,l,....,z-/-l,z-/ + l, 

J2 Pk0R)n X Pi{-dR)T,, 

k+l=n 
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(5.32) 


PK5«)[0,1,2,...,2- 


Mr = 


[0,1,2,..., z 


-2,z + /- 


(5.33) 


FK-^i?)[0,l,2,...,2- 


Mr = 


(-i)'[o,i,.. 


.,z - / - 1, 


where 














C/co,0 


Cfco,l 




f5 34^1 


kn k- ,]r — 


Cfci.o 




Cfci,i-1 








• • • 








Co,fco 


Com 




(5.35) 


[ko, h, /C2, . . . , fci-i]/? = 


Ci,fco 


Ci,fci 


Ci,fc,_i 








Cj-i,fci • • • 


Cj-i,fci_i 



where Pk are Schur polynomial as defined in f l2.10p . 

In fact. Lemma [5.41 is a special case of abstract general formula of Schur function [18] 

(5.36) SY{d)T^ = TY{t). 

Here := [0, 1, 2, . . . , i — 1] is the standard Wronskian determinant, Y := (Yq, Yi,Y2, . . . , Yi^i) 
and 



S- 



Y 



P 



Yi- 



Y^-2-l 



Pyi+1 Pyo 

Pyi Pyo-1 

PYi-i+3 PYa-i+2 

PYi-i+2 Pyo-i+I 



Y = (Yq, Yi,Y2, . . . , Yi_i) (Yq > Yi > Y2 > . . . > Yi_i) corresponds to Young diagram with Yq 
boxes at the first row, Yi boxes at the second row and so on. Notice formula [18] 



(5.37) 



Sy{-t) = (-l)l^l5y,(t), 



where Y' is the conjugate Young diagram of Y. Eq. fl5.36p together with eg. (15.371) further leads 
to the lemma above easily. 

Then we define wave functions Wl = {Wli, Wl2, ■ ■ ■), = {Wri, Wii2, ■ ■ -Y and Wr = 
{WRi,WR2,...f with 



(5.38) WLiiX^,t) 



Ti-l 



Co,o C'0,1 
Ci,o Ci,i 



Co,i-2 1 
C14-2 Aiv 



(5.39) WR,{X^,t) 





Co,o 
Ci,o 


Co,i 

Ci,i . 


■ ^0,j-l 

■ Cij-i 




Cj-2,0 
1 


Cj-2,1 ■ 


■ Cj-2,j-l 
A M 
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Co,o 


Co,i 


• Coj-i 




Ci,o 




■ Cij^i 






Cj-2,1 . 






1 


Am 


A M 



(5.40) WR,{\M^t) 



which satisfy the following orthogonality relation, 

(5.41) {WLi, Wr,) = {Wu, Wr^) = 

where 

and inner product (, ) of functions A and B is defined as 

{A,B) := j j p{X,^)S{X^ - fi^)A{X,t)B{fi,t)dXdfi. 

Therefore tau functions have another form as 

Tm := (let ({Wu,WRj) 



l<i,j<r7i 

The entries of the matrix representation of the Lax operator C can be then calculated by 



Pi-j+NiDljTjTi^l 



Ti-lTj 



[0,1,. 


. j — N + m 


-2,3- 


N + m 


Co,o 


Co,i 


Co,i-2 


1 


Ci,o 


Ci,i . . . 




Xn 




. . . 




X N 



m+n=i—j+N 
^ i-j+N 

V [0,l,...,j-2,j + m-l]^ 

* ^ 3 m=0 



So 

(5.42) a,,, = (Al^i„W^^,) 

which are given by the matrix representations of the eigenvalue problems CWl = XWl and 
WrC = XWR^see [12j for the details). Till now, we have solved the BTH using orthogonal 
polynomials. 

If we denote VVl and VV^^ as corresponding matrix forms of Wl (15.151) and (15.161) 
respectively, then Wl and Wr can be represented by matrices Wl and VV^^ respectively as 
following. 
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Co,o 


Co,i 


• Coj_i 


eCi?.{A^,t) 


Ci,o 


Ci,i . 


• Cij_i 




Cj-2,0 


Ci-2,1 • 


• Cj_2,j-1 




1 


Am 


A M 



Because 



A 



we can get 



Wj. 



( \ \ 

Xn 

V ■ y 



/ 1 \ 



V ■ / 



/ 1^ \ 

Aiv 

Ajv 

V ■ J 



( 1 \ 



Pi(x,t,A) 



Aiv 
Aiv 



V ■ / 



1 



T2 



T2 



\ I 

where VL{x,t,K) is as matrix (13.361) . This agree with the definition of Wli in eq.f l5.38p . Also 
similarly we can get 



R 



( 1 All 



Am 



^tMlL + Ai^^ 



(1 A 



Am 



P2{-dR)T2 , Pl(-8fl)r2 
T3 -ra 



Am +21 A 



T3 



where A) is as matrix fl3.39p . This also agrees with the definition of Wm in eq. (l5.39p . 

Formal factorization mentioned above is about infinite-sized Lax matrix. In the next section, 
we will consider its finite-sized truncation. Then we can find finite-sized Lax matrix is in fact 
nilpotent because of dressing structure fl3.21l) . This finite-sized case corresponds to rational 
solutions of the BTH which will be considered in the next section. 

6. Rational solutions of the (A^, M)-BTH 

It is well known that the r-function of the original tridiagonal Toda lattice, i.e. (1, 1)-BTH, 
has the Schur polynomial solutions associated with rectangular Young diagrams. So what kind 
of Young diagrams correspond to the BTH become an interesting question. In this section, 
we only consider homogeneous rational solution of (A^, M)-BTH which is one kind of most 
interesting solutions in nonlinear integrable systems. 

In order to describe homogeneous rational solution of (A^, M)-BTH (A^ < M), we firstly set 
deg{tifi) = degitofi) = MN. Then we get deg{Pm{ta)) = inriM and deg^Pnit/^)) = nN . 
Define 



m,n 



and choose the following special homogeneous polynomials Pkita^tis) with degree k as ri, 

N-l M-l 

Pk{ta,tp) ■=^^ c'^,q Pm{ta)Pn{tl3), 



p=0 q=0 mM+nN=k 
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where ta and denote 

ta = {ta,n ■■ I < a < N,n = 0,1,2, .. .} 

tp = {W : -M + 1 < /3 < 0, 72 = 0, 1, 2, . . .}, 

and k can be any number in the set {niM + nN;m,n G Z+}. The polynomials Pm{ta) and 
Pnitp) are the elementary Schur polynomials, and they satisfy the following relations, 



dPM 



P 



Q, - m-Af(m' + l)+a'-lvay) 

Note here that {c^^^lO < m,p < N — 1,0 < n,q < M — 1} can be arbitrary constants. 



dPn{tp) 



P 



n-M{n'+l)-l3' 



Define 



Af-l M-1 



p: 



LV 



EE E 



c: 



m+l,n+l' 
P,<1 



Pm{ta)Pn{ti3), 



p=0 q=Q mM+nN=k 

and the rational solutions for (N,M)-BTH have the following diagram representation, 

Dj = {k- {j - 1)N, k- {j -2)N - M,...,k- {j - 1)M}, k = 0,N, M, 2N, N + M,2M,.... 

The difference between two adjacent two numbers in Dj is M — N. 

We denote the tau function corresponding to Dj{k) as Tj^Oj which have the following form 



'Tj,D,(k) - So^ik) 



P, 



P 



(0,0) 

5: 

(0,1) 

k-N 



p{1.0) 
^k-M 

p(l>l) 

^k-M-N 



^k-(j-l)M 



p 



0-1,1) 



p 



(Oj-l) 



p 



{i,i-i) 



k-(j-l)N k-M-{j-l)N 



P 



k~(j~l)M~N 

0-ij-i) 



k-{j-l){M+N) 



5(0,0) 



where t,- denotes the j-th tau function(j x j determinant) generated by -P^ • The range 
of rank j depends on the choice of k. This kind of diagram like Dj is not classical Young 
diagram. It is a kind of generalized diagram which counts the homogeneous degree which comes 
from the multiplication of two classical Shur functions. We can call this kind of generalized 
diagram degree diagram. Same as Young diagram, the tau functions represented by degree 
diagram is also Wronskian form, the derivative is about dt_j^f^j^Q or dt^^g. Because the scale of 
degree in definition before(e.g. deg{Pm{ta)) = rnM) is bigger than common degree of Schur 
polynomial((ie5f(Pm(t)) = m), the difference of subscript between two adjacent rows(columns) 
is N{M) not 1(1). From this point, it is also different from Hankel determinant. 

In the following we firstly only consider the case when and M are co-prime integers. When 
they are not co-prime, just divide them by the GCD of them and use the theory of co-prime case 
in the following. In fact, we can find for fixed value oi p, q, all the coefficients of a homogeneous 
polynomial will be the same because they all satisfy relation f l5.25p . Then we get 

N-l M-1 



EE E 



c. 



m+l,n+l' 
PS 



Pm{ta)Pn{tl3) 



p=0 q=0 mM+nN=k 
N-l M-1 

^p^q Prriiita) Pfii^tj^) -\- Cp^ Pm2it oi) Pn2it p) 

p=0 g=0 



EE 



•P,9 
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fni+l,ni+l' 



PmitajPnitls) 



p=0 5=0 



mM+nN=k 



mM+nN=k 



where 



c, 



,mi+i,ni+/' 
'P>9 



= C, 



p,q 

N-l M-1 

EE' 



,mi+(,ni+r 



niiM + rijA^ = k, m. 



p=0 5=0 

nii - N, nj+i 



M, l<i< l{k), 



and number l{k) depends on k. For simplicity, here we just consider the case with all coefficients 
equal 1 which is our central consideration in this section and define 

Pk = Pm{ta)Pn{tl3) = -Pmi (^a)-fni (^^) + Pm2{ta)Pn2{'tp) + " 1" -Pm((fc) (^a)-Pn((fc) (^Z?) ■ 

mM+nN=k 



Then the r-function generated by Ti 
determinant, 



Pfc can be expressed by a double- Wronskian 



Ts{k,ta,ti3) 



Pk 
Pk-N 



_Pk-M 
Pk-M-N 



Jk-{s-l)M 
Pk-{s-l)M- 



N 



Pi 



k-(s-l)N ^k-M-(s-l)N 



^ Pnii (tc 



P. 
P 



mi — 1 



mi— 2 



Pm2 (^a) 
Pm2~l (j^ce) 
Pm2—2ita) 



k-is-l){M+N) 



mj — 1 (^a / 
^m(-2( 



Pm,-2{ta) 



Pmi—s+l{^a) j 



X 



y -fmi— s+l(^a) -fm2— s+l(^a) 

/ Pniitp) Pni-litp) Pni-2{tl3) ■■■ -Pni -s+1 (^Z?) \ 
Pn2{'tl3) -Pn2-l(^^) Pn2-2{tl3) ■■■ Pn2-s+l{t0) 



\ Pniitis) Pni-litp) Pni-2{tl3) 



ni—s+1 



Conversely, for a fixed size j of Lax matrix for [N, M)-BTH, the choices of k for ti is in set 

(6.1) Kj := {k\k = {j - 1)NM + mM + nN,m,n e Z+,0 < m < N,0 < n < M} . 

We can see that the number of elements in set Kj is A^M(When they are not co-prime, this 
number will be j^^, where (iV, M) is CCD of N and M). When the values of N, M, j, m, n 
are chosen, a series of non-vanishing tau functions corresponding to them will be fixed. 

In the following, we will consider rational solutions of (A^, M)-BTH with finite-sized Lax 
matrix. For (A^, M)-BTH, the size of minimal Lax matrix is (M + 1) x (M + 1). This minimal 
Lax matrix has M non-vanishing r-functions and anyone's degree has M — N jumps between 
adjacent rows. For the special case oi N — M, degree diagrams are always rectangle which can 
also be seen from definition of Dj. 
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Besides considering the degree diagrams, it is more interesting to consider the decomposition 
of degree diagrams into representation of Young diagrams. In fact the Young diagram represen- 
tation of general BTH has a form of muhiphcation of two different groups of Young diagrams 
which will be shown in the following theorem. 

Theorem 6.1. For j x j{j > M + l)-sized Lax matrix of {N, M)-BTH (denoted as {N, M)j,,j), 
after choosing the value of k as {j — 1)MN + mM + nN, Young diagram representation of a 
series of corresponding tau functions are as following 

Tl = 'S'{(j-l-a)Af+m)(^a)5'(„+aM)(^/3)) 
0<a<j-l 

it,), 

0<a<b<j-l 

T3 = S((j-i-a)N+m-2,(j-l-b)N+m-l,(j-l-c)N+m) it,), 

0<a<b<c<j-l 

'^s = ^ ^ S(^(^j^i^ai)N+m-s+l,{j-l-a2)N+m-s+2,...,(j-l-as)N+m){ta) 

0<ai<a2<---<as<j — l 

S[n+asM —s+l,...,n+a2M—l,n+aiM) it,), 
'^j ~ S[(j^i)(^N-l)+m,{j-2)iN~l)+m,...,m)ita)S(n+{j-l)(M-l),n+{j-2)iM~l),...,n)itl3)- 

Proof. To prove this theorem, one need to use Cauchy-Binet formula. The process is quite 
complicated because of huge sizes of matrices. So we will omit the proof. One can understand 
the patten by the following example, i.e. (2, 3)-BTH. □ 

To see it clearly, we give some specific examples in the following. 

Example 6.2. {2,3)-BTH, has 6 sets of t -functions for each size of Lax matrix and the degree 
diagram for every tau function has one jump between adjacent rows. See (2,3)4x4 i'^ detail as 
following degree diagram 



(6.2) (2,3)4x4 



'(0,0) 


^l,{18} 




T2,{16,15} 


-)■ 


''"3, {14,13,12} 




'r4,{12,ll,10,9}5 


(0,1) 


Tl,{20} 


— 


'i^2,{18,17} 


-> 


T3, {16,15, 14} 


— 


'i^4,{14,13,12,ll} 


(1,0) 


ri,{21} 


— 


T2,{19,18} 


— )■ 


'^'3, {17,16,15} 


— )■ 


''4,{15, 14,13,12} 


(0,2) 


■'"l,{22} 




''"2, {20,19} 


— > 


''"3,{18, 17,16} 


— > 


''"4,{16, 15, 14,13} 


(1,1) 


^,{23} 




'''2, {21, 20} 


— > 


7"3,{19,18,17} 


— > 


''"4,{17,16,15,14} 


.(1,2) 


n,{25} 




'''2, {23,22} 




''3,{21, 20,19} 




'74,{19, 18, 17,16} 



where {ip,q),0 < p < 2,0 < q < 3} denote the value of {m,n) in value of k, i.e. (16. ip . Here 
k takes values in {18,20,21,22,23,25}, i.e. the values in bracket o/ri^{}. Every tau function 
generates a series of tau functions which are connected by right arrow in (16.21) . r^^j ^...^ } 
represents the l-th tau function whose degree diagram is in the bracket {,..., }. In (16.21) . 
we can use product of Young diagrams to represent the four tau functions of the first line, i.e. 
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(0,0) case as following by Cauchy-Binet formula 

( 1 \ 

S I I I I I I I {ta)Sti){tp) + S rrm {ta)S ojj {tp) + S an (tn)S rm-m 



Tl,{18} 



( Peita) Piita) P2{ta) 1 ) 



{tp)+S^{ta)S 







T2 = T2,{16,15} = 



P(^{ta) PA{ta) P2{ta) 1 
P^{ta) Psita) Pl{ta) 



/I \ 

S {ta)S U2 {tp) + S gxn {ta)S rrrm (^/s) + 5 mTTI (ia)'S' (i^) 



/ 1 \ 

P3(t/3) P2{tp) Pl{tp) 

Peitp) P^itp) P^itp) 

\ P^m) Ps{t^) P,{tp) ) 



TS — T3,{14,13,12} — 



Peita) P^ta) P2{Q 1 
PM Pzito) Pl(ta) 
P4{ta) P2{ta) 1 



T4 — r4_{i2,ll,lo,9} — 



/ Peita) PM P2{tc 

P^iQ PsiQ Pl{tc 

PA{ta) P2{ta) 1 

\ PsiQ P,{Q 
5' pi (t«)5' pFp (tff). 



1 \ / 1 
















\ 

Psitf,) P2{tp) Piit^) 1 
Peitp) P^itp) P^itp) P^itp) 
\ P,{tp) Ps{t^) Pj{^) P,{tp) I 



After general theory on homogeneous rational solutions of the (A^, A'/)-BTH, as a special but 
important case, rational solutions of the (1, M)-BTH will be considered in the next subsection. 

6.1. Rational solutions of the (1,M)-BTH. Since the ti^n flows are same as the to,n flows, 
we use tl „ + to,n as a new variable and identify as to,n- Then the rational solutions for the 
(1, M)-BTH are obtained from the r-function. 



(6.3) T^{k,t) 



Pk 

Pk-l 



Pk~M 
Pk-M-l 



-< fc-(j-l)M 
-Pfc-(j-l)M-l 



fe-0-l)(M+l) 



where P„ = Pnitp) with to,n = ^o,n + ^i,n- This r-function can be given by the Schur polynomial 
associated with the Young diagram which is same as degree diagram mentioned above for 
(1,M)-BTH, i.e. Tj{k) = Sy^ik) with 



Y^{k) = {k-3 + l,k-3 + l-{M-l),...,k-{j-l)M) for j = l,2,...,l + 
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k_ 
M 



where [-^J denotes the biggest integer which is less than or equal Note here that the number 
of boxes in the diagram increases by M — 1 between adjacent rows. Let us now give some 
examples of the r- functions for the (1, M)-BTH with specific size r of the Lax matrix, denoted 
by (1, M)rxr- For a given size r(> M), the choices of k are in the set {{r — 1)M, . . . , rM — 1}, 
that is, there are M choices of the first member of the r- functions, ti = Pk- Then each value 
of k generates r tau functions ordered from ti to r^. 

Example 6.3. {1,1)-BTH, i.e. the original Toda lattice, has only one set of t- functions for 
each size of Lax matrix; 

(1, 1)2x2 : Tn T^, 
(1, 1)3x3 : Tm rg r^, 
(1,1)4x4 : r nn rg r^, 



(6.4) 



where = 1 . 

{1,2)-BTH has two sets of t -functions, and the Young diagram for each r-function has one 
jump between adjacent rows: 



(6.5) 



f (1,2)3x3 

1, 2)4x4 



r '^ffffP ~^ 

T — )■ r 1 1 1 1 1 1 1 T 



□ -)■ r 



Similarly {1,3)-BTH has three sets of tau functions, and the Young diagram has two jumps 
between adjacent rows. 



r(i,3) 



4x4 



: T — >■ r | 1 1 1 1 1 1 1 1 — >■ T 

r — >■ T | 1 1 1 1 1 1 1 1 1 — >■ r 

T -n — >■ TI I I I I I I I I □ — >■ T 



n — >■ T 



m — >■ r 



m — )■ r 



7. Conclusions and discussions 

We proved the equivalence between (A^, M)-BTH and (M, A^)-BTH, derived the primary Hi- 
rota equations of the (A^, M)-BTH, and found several explicit formulas about solutions for the 
BTH using orthogonal polynomials in the matrix form. We also constructed some rational solu- 
tions of the BTH which are parameterized by the products of Schur polynomials corresponding 
to non-rectangular Young diagrams. It may be interesting to find their significance in terms of 
the representation theory, as in the case of the original Toda lattice where the rational solutions 
are given by the Schur polynomials of rectangular Young diagrams and they are the Virasoro 
singular vectors. 
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